Magnetic field-induced control of transport 
in multiterminal focusing quantum billiards 
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By exploring the four-terminal transmission of a semi-elliptic open quantum billiard in depen- 
dence of its geometry and an applied magnetic field, it is shown that a controllable switching of 
currents between the four terminals can be obtained. Depending on the eccentricity of the semi- 
ellipse and the width and placement of the leads, high transmittivity at zero magnetic field is reached 
either through states guided along the curved boundary or focused onto the straight boundary of 
the billiard. For small eccentricity, attachment of leads at the ellipse foci can yield optimized corre- 
sponding transmission, while departures from this behavior demonstrate the inapplicability of solely 
classical considerations in the deep quantum regime. The geometrically determined transmission is 
altered by the phase-modulating and deflecting effect of the magnetic field, which switches the pairs 
of leads connected by high transmittivity. It is shown that the elliptic boundary is responsible for 
these very special transport properties. At higher field strengths edge states form and the multiter- 
minal transmission coefficients are determined by the topology of the billiard. The combination of 
magnetotransport with geometrically optimized transmission behavior leads to an efficient control 
of the current through the multiterminal structure. 
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I. INTRODUCTION 

Open particle billiards are two-dimensional (2D) 
scattering structures that confine ballistically moving 
(charged) particles within a region of space of certain 
geometry, at the boundary of which openings allow 
the particle to escape into leads. Their quantum ver- 
sion models nanoscale transport devices, usually real- 
ized by controllably confining a 2D electron gas at a 
semiconductor interface, and sets the grounds for the 
theoretical description and experimental investigation 
of coherent transport in the mesoscopic regime. The 
transport properties of such billiards are extensively ex- 
amined in the context of various interesting phenom- 
ena, including the quantization of conductance in meso- 
scopic systemSfii^ the appearance of Fano resonances 
in transmission;^"— resonance trapping^i^ shot noise in 
quantum dots^ii^ Andreev tunneling and reflection ;^^!^^ 
as well as conductance fluctuations r^iii^ localization 
effectsr^"— decoherenc o^'^^1^^ and dephasing^i^S in bal- 
listic nanostructures. Billiards also serve as ideal sys- 
tems for the study of the relation of quantum transport 
to its classical counterpart and the crossover between 
them^i^^-^ 

In all these phenomena the quantum nature of trans- 
port manifests itself in interference effects, which arise 
due to the phase coherence of the underlying scatter- 
ing process. The interference of the scattering states 
and its impact on the inter-lead transmission of a quan- 
tum billiard is determined by its geometry including 
the placement of the leads on its boundaryi^ For the 
(semi-) elliptic geometry considered here, the classical 
dynamics of the closed billiard is regular, with ballis- 
tic particle trajectories divided into so called librators 
and rotators, which intersect its major axis at the seg- 
ment between the foci and the segments between the 
foci and the boundary, respectivelyi^ Librators and rota- 
tors correspond to quantum eigenmodes^^ localized about 



the minor semi-axis (also called 'bouncing ball modes') 
and along the elliptic boundary (also called 'whispering 
gallery modes'), respectivelyi^ Attaching leads to the 
straight boundary of the semi-ellipse results in a gener- 
alized open mushroom^! (Bunimovich) billiard with mul- 
tiple stems of infinite length. The chaotic trajectories of 
the closed mushroom, entering its stems, escape into the 
leads in the open billiard and contribute to transmission. 

In the presence of a magnetic vector potential the 
Aharonov-Bohm (AB) effeclj^ induces intricate changes 
with respect to the interference patterns. Generalized 
from its occurrence in quantum ring a^^'^° to spatially ex- 
tended billiards^^^iS it leads to fluctuations of the trans- 
mission coefficients with the magnetic field strength. At 
high field strengths, where unconfined electrons would oc- 
cupy localized Landau states, in the billiard edge states 
form24 which mediate the propagation along the bound- 
ary, with their interference at the lead openings leading 
to oscillatory transmission behavior.— 

In the present article we investigate the transmission 
behavior of a 4-terminal semi-elliptic quantum billiard in 
dependence of its geometrical characteristics and examine 
how an injected wave can be guided selectively to a chosen 
output terminal using a magnetic field. It is shown that 
the zero-field transmittivity between leads attached to the 
billiard highly depends on the accessibility, through their 
coupling to these leads, of librator and rotator modes. 
The crucial role of the convex billiard boundary in the 
transmittivity of such states has been described in Ref. [2^ 
for a 2-terminal semi-circular billiard. We in turn calcu- 
late the 4-terminal transmission coefficients for varying 
eccentricity of the semi-ellipse, exploring the dependence 
of the multilead transmittivity on the curvature of the 
boundary. The crossover from librator to rotator modes 
being the dominant transmission mediators is revealed by 
altering the lead positions along the straight part of the 
billiard boundary. It is shown that an optimal transmit- 
tivity between both pairs of symmetrically placed leads 
can be achieved by an appropriate choice of eccentricity 
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FIG. 1. (Color online) Geometry of the open billiard with 
indicated length parameters and lead labeling for eccentricity 
e — f/a = \J\^^W-la?' — 0.35, lead width w — a/5, and inner 
lead positions ±Z = ±a/2. The two dots at {x, y) = (±/, 0) are 
the foci of the (semi-) ellipse and the dashed lines correspond 
to alternative setups (see text). 



and lead positioning. For small eccentricity, the place- 
ment of attached leads at the foci of the ellipse yields 
a high corresponding overall transmission coefficient, as 
would be expected classically. As a result of interference, 
though, this condition docs not apply for a generic setup, 
which implies a departure from purely classical considera- 
tions. The importance of the rotator and librator modes 
is further assessed by gradually perturbing the billiard 
with a circular disk placed on the curved boundary or in 
the interior. 

With restrictions due to symmetry, the magnetic field 
changes each transmission coefficient differently; the 
question thus arises whether the multiterminal setup can 
function as a selective switch between the terminals by 
tuning the magnetic field, in the sense that transmission 
of an incoming particle is efffciently favored to specific 
leads and suppressed to others. Calculating the multi- 
terminal transmission coefficients of selected setups for 
varying magnetic field, we show that such output con- 
trollability is indeed achieved: it results from the form of 
the scattering wave function determined by the billiard 
geometry and the magnetic fieldj^ 

The paper is organized as follows. In Section |lT] the 
geometrical setup of the 2D billiard is specified, the theo- 
retical framework for quantum transport is outlined and 
the computational approach is presented. Section IIII Al 
summarizes the consequences of the symmetries of the 
system with respect to the multiterminal transmission. In 
Section UlI BI the main features of the obtained multilcad 
transmission spectra are discussed, along with a descrip- 
tion of the underlying mechanisms. This is followed by an 
analysis of the mean transmission components in depen- 
dence of the geometric properties of the billiard in Section 
IIII CI The impact of the magnetic field on transport is 
discussed in Scction flll Dl concluding on the induced con- 
trollability of transmission to selected output leads. The 
necessity of the billiard properties for controllable com- 
bined output is demonstrated in Section [ill El and finally 
its functionality is shown in Section IIII Fl in terms of the 
linear conductance at finite temperature. Section |IV] con- 
tains a brief summary of our conclusions. 



II. SETUP AND COMPUTATIONAL 
APPROACH 

The geometry of the 4-terminal hard-wall billiard is 
show n in Fig.[Tl It consists of a semi-ellipse of eccentricity 
e = \/l — \P- ja? , where a and h are the major and minor 
semi-axes respectively, on the straight border of which 
four vertical semi- infinite leads of equal width w < a/2 
arc attached, symmetrically about the minor axis of the 
ellipse a; = 0. The elliptic boundary is smoothly contin- 
ued into the outer leads 1 and 2, between which the inner 
leads 3 and 4 are centered at distance I from the origin. 
Semi-infinite leads attached to the closed billiard repre- 
sent the coupling to particle reservoirs, from which there 
is no reflection back into the billiard. 

To calculate the multi-terminal transmission coeffi- 
cients of the system, Dirichlet boundary conditions are 
imposed on the scattering wave function along the bound- 
ary of the billiard, defined by a hard- wall potential 
V{y)^ Adapting to natural units, we set h = —e = 
m = 1, where e < and m denote the charge and effec- 
tive mass, respectively, of a spinlcss particle. The single- 
particle Hamiltonian is then written 



H 



-(k 

2 ^ 



■AY + V 



(1) 



where the vector potential A generates a magnetic field 
B = V X A = Bz perpendicular to the plane of the 
structure, with strength B that is homogeneous over the 
extent of the billiard and falls off linearly to zero over 
a length w 30w in the attached leads. The particle is 
incident in one of the four leads with energy 
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where n = 1, 2, ... labels the subbands of the longitudinal 
momentum ky along the unperturbed leads, generated 
by the transversal confinement to their common width 
w. The scaled momentum k = k ■ w/n ~ \J2E ■ w/t: 
thus varies continuously in the interval n < k < 7i + 1 for 
motion in the n-th subband. 

The Hamiltonian is discretized on a square grid of 
unit lattice constant in the tight-binding approach, with 
the magnetic vector potential incorporated through the 
Pcicrls substitutionr^ Considering (spin polarized) elec- 
tronic transport at a GaAs/AlGaAs (m = 0.069 nif.) in- 
terface and setting the lattice constant unit to a = 2 nm, 
the unit of energy becomes h/mo? = 0.276 meV and 
the unit of field strength h/\e\a'^ = 164.55 T, while the 
lengths in the system are scaled by a reference ellipse 
major semi-axis a = 100a = 200 nm. 

The coupling of the system to the external semi-infinite 
leads i = 1, 2, 3, 4 attached to the system is described by 
self-energies S^, which are analytically obtained for B = 
and contribute non-Hermitian blocks to the Hamilto- 
nian matrixi^ From the single-particle (retarded) Green 
function of the system 



where 



G(^) = [M- (H + S)]-i 
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the part describing the propagation from lead j to 
lead i is computed using a parallel implementation of the 
recursive Green function method (RGM) , where a decom- 
position scheme among communicating processors allows 
for the computation to be done in a parallel mannerj^ 

Since in general B ^ 0, the self-energies are addressed 
to grid points in the leads where the field strength has 
decreased to zero. The computation of G is then more 
efficient if the system is decomposed into parts, so called 
modules, of increased symmetry, by utilizing our own im- 
plementation of a modular version of the RGM, intro- 
duced in Ref. H^- The considered setup is thus assembled 
by two types of modules: the semi-elliptic scatterer at 
constant magnetic field strength B, and a lead part with 
the magnetic field decreasing linearly to zero. The mod- 
ules are then subsequently connected to form the com- 
plete scatterer, and the Green function of the connected 
modules is obtained in each step by solving a matrix 
Dyson equation in tight-binding form4£ 

The scattering matrix elements connecting the leads 
of the system are expressed in terms of the Green func- 
tion by the Fisher-Lee relation;^ which leads to a com- 
pact trace formula for the multiterminal transmission 
coefficients 
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TABLE L Multiterminal transmission coefficients Tij{B) from 
lead j (columns) to lead i (rows) , deduced from the coefficients 
721, 132, Ti3 and T-m at magnetic field strength B (see text). 
The surrounding boxes of the Ti^tjj correspond to the plotted 
line types in Figs. |3]-|H1 



Time reversal symmetry (TRS) yields a Hermitian conju- 
gated scattering matrix under inversion of the magnetic 
field, implying for the transmission coefficients)^ 



T,,{E;B)^Tj,{E;-B) 



(8) 



r,,(ii;) = Tr (r.Gr.Gt) (t^j) 



(5) 



with Tj = — Sj). For i ^ j, the reflection coefficient 
of each lead j, Tjj = Rj, is given by the sum rule 

n{E) =Y,T,,{E) =Y,T,,{E) 

i i 

= J2t^,{e) + r,{e) , (6) 

where 'n(E) ^ iTit[K,{E)] is the number of open channels 
in the leads (of common width w) at energy E, int[ ] de- 
noting the integer part. This sum rule results, for coher- 
ent transport, from the unitarity of the scattering matrix, 
which in turn is a consequence of probability flux conser- 
vation. 

The Green function of the system is also used to calcu- 
late the local density of states (LDOS) at site r through 
the relation 



p(r,ii;) = ^(r|*(ii;)|r) 



(7) 



where $ = GFG^ is the spectral function and F generally 
a weighted sum over the F^ according to the distributions 
of incoming states in the leads (the Fermi distribution in 
the case of electrons). In the cases presented here, we 
have chosen F = F^ for a certain lead i, so that p{r,E) 
corresponds to the probability density resulting from an 
incoming monochromatic wave of energy E in lead i. 



This reciprocity relation halves the number of indepen- 
dent transmission coefficients Ti^j^j when both field di- 
rections are considered. The rcfiection coefficients Rj = 
Ti=jj are given by the sum rule Eq. ^ and remain the 
same under field reversal [cf. Eq. ([8|)], reducing the num- 
ber of independent coefficients by the number of leads. 
For a 4-terminal billiard the 4 x 4 = 16 coefficients are 
thus reduced to 6 independent ones. 

Iw our billiard the spatial reflection symmetry about 
the j/-axis introduces additional relations between sym- 
metric pairs of leads. If leads i (at Xi) and j (at Xj) are 
placed symmetrically to leads i' (at —Xi) and / (at — Xj), 
then 



T,,{E;B)=T,,j,{E;-B) 



(9) 



as the equations of motion for a (charged) particle are 
invariant under the transformation (x, B) — > {—x, —B) 
in the symmetric billiard. Explicitly, we get the two 
additional relations T24iB) = Ti3{-B) and T4i{B) = 
T'i2{—B), reducing the number of independent coefficients 
to 4 (if i' = j and j' = i, this reflection symmetry coin- 
cides with the TRS). 

In the following, we will work explicitly with the coef- 
flcients T21, T32, T13 and T^i with input in each one of 
the four leads, because this set serves best for our discus- 
sion of the results. In Table U all transmission coefficients 
Tij{B) are explicitly expressed in terms of these four. 



B. Transmission spectra at zero magnetic field 



III. RESULTS AND DISCUSSION 
A. Symmetries of the transmission coefficients 

Before investigating the multiterminal transmission co- 
efficients in varying geometry and field, we show how 
symmetries present in the system can be used to reduce 
the number of independent coefficients to be calculated. 



The generic features of the zero-field multiterminal 
transmission spectra are presented for a geometric setup 
with relatively high overall transmission between inner 
leads, since this will prove to be a key property for the 
desired controllability of output terminal. The (k) are 
shown in Fig.[2JA) for k within the first channel (n = 1). 
All terminal combinations are represented, since their 
relations in Table |T] simplify accordingly for B = —B. 



4 



1 1.2 1.4 1.6 1.8 2 
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FIG. 2. (Color online) (A) The four independent transmission 
coefficients Tij at zero magnetic field, as a function of the 
scaled energy k = k ■ w/n = \/2E ■ w/tv, for the elliptic (e — 
£3 = 0.35, black solid line) and circular (blue dotted line) 
billiard border, both with w — a/5 and I = esa. (B) LDOS 
for the semi-elliptic billiard at different energies (a), (b), (c), 
(d) indicated by the vertical lines in (A), with the particle 
incident in the (i) outer and (ii) inner left lead. The colormap 
scales with ^y p{x, y) from white [p — 0) to black (p — max) 
and is normalized to its maximal value in each plot. 



Fig. m^B) shows the LDOS at selected energies for dif- 
ferent leads of incident wave. Note that the presence of 
interference fringes in a lead signifies back-reflection into 
that lead. 

As we see in T2i(k), transmission is overall close to 
unity between the outer leads. This results from rotator 
modes of the semi-elliptic billiard that would leak into 
finite stems at the outer lead positions, and therefore are 
strongly coupled to these leads in the open system. These 
leaking states constitute non-resonant pathways for trans- 
port, whose superposition leads to a high transmission 
background in T21 [see LDOS in Fig.mB)(i,d)], smoothly 
varying m energy] 1 

States weakly coupled to the outer leads [see LDOS 
in Fig. [2lB) (i,a) and (i,b)] constitute resonant path- 
ways, whose interference with the non-resonant pathways 
leads to sharp resonances in T2i{k) of width proportional 
to their coupling, which possess the characteristic Fano 
lincshape asymmetry4^"— In the case of a single non- 
resonant pathway, the asymmetry is caused by a trans- 
mission zero close to the resonant energy, owing to com- 
plete destructive interference between the resonant and 
non-resonant state. In our system each quasi-bound state 
(resonant pathway) in general couples to multiple multi- 
terminal leaking states (non-resonant pathways), which 



renders the total interference partially destructive and 
thus raises the minimum of the Fano resonance from 
zero4i Corresponding to the eigenstates of the closed 
semi-ellipse, the resonances superimposed on the trans- 
mission background appear in series of different quasi- 
periodicity in k, determined by the quantization of the 
wave- number of the semi-elliptic modes (in analogy with 
the detailed description in Rcf. |42 for the oval billiard). 

The coefficient T34 is also overall high in the first chan- 
nel for the chosen eccentricity and lead positioning, in 
this case resulting from the strong coupling mostly of li- 
brator modes to the inner leads. The convexity of the 
boundary plays a crucial role for this behavior, since it 
focuses the scattering wave function, incident in an inner 
lead, around the middle of the straight boundary of the 
biUiard [see Fig.[2jB)(ii,c) at which energy T34 practically 
reaches unity]. Note that, since the leads in this setup are 
centered at the foci of the ellipse, classically both libra- 
tors and rotators intersect the lead openings; quantum 
mechanically though, there is a larger number of eigen- 
modes of librator type with maxima at the focii^ Due to 
interference between these modes, setting Z = / = ea is 
not a necessity for the acquired high inner-lead transmis- 
sion; it depends also on the chosen e, as will be shown in 

Section inrci 

In the sharp resonant dips in the high background 
arc at the same positions as in T21 but generally of dif- 
ferent width, arising from the same quasi-bound states 
coupling with different strength to the inner leads. Also, 
the Fano minima in T34 are of different height than in 
T21, since different non-resonant transport paths are pro- 
vided by the librator modes, interfering with the reso- 
nant states. For some sharp resonances in T34 (e.g. at 
K sa 1.145 and 1.365) the Fano asymmetry is more dis- 
tinct, as they lie within a dip of the transmission back- 
ground. 

As a consequence of probability flux conservation 
[Eq. ([5])], unit transmission in either T21 = T12 or T34 = 
r43 leads to vanishing transmission in both r32 = T23 = 
T41 = Ti4 and T13 = T^i ~ T24 = T42- Thus, transmis- 
sion between an outer and an inner lead, represented by 
the coefhcicnts T32 and ri3 in Fig. A) , is almost zero 
over the whole channel for the chosen geometric parame- 
ters. It exhibits resonant peaks, coinciding with dips for 
the symmetric lead pairs, albeit of rather low amplitude, 
since the the incoming wave is mostly reflected into the 
same lead [see Fig. [5][B)(a)] or transmitted to the sym- 
metrically placed lead [see Fig. [2jB)(b)]. As each eigen- 
state of the semi-elliptic billiard is symmetric, it gener- 
ically couples with different strength to inner and outer 
leads; desymmetrized lead positions then lead to lower 
transmittivity of the corresponding resonant states in the 
open systemi^ 

When the eccentricity of the billiard is slightly changed, 
the eigenstate wavelength in the elliptic coordinates is ac- 
cordingly modified^^ and consequently the corresponding 
transmission resonances shift in k. This is evident in 
Fig.[2)JA) for the semi-circular billiard (dotted cyan line): 
its area is slightly larger than the semi-ellipse, so that 
the resonances are shifted to lower k. The dips in r2i 
and r34 (or peaks in T32 and T13) are overall broader for 
the semi-circular billiard, owing to the enhanced coupling 
of its eigenstates to both inner and outer leads. It also 
shows a pronounced imbalance between r32 and T13 at 
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FIG. 3. (Color online) Zero-field channel-averaged transmis- 
sion coefficients Tij between terminals = (2,1) (red 
dashed-dotted line), (3,4) (blue solid line), (1,3) (cyan dot- 
ted line), (3, 2) (green dashed line) with lead widths and po- 
sitions as in Fig. [2] and eccentricities (from top to bottom 
row) ei = 0.78, £2 = 0.65, £3 = 0.35, for the semi-elliptic 
billiard with a disk inserted in the billiard (a) centered at 
(2;, y) = (0, b) with radius r varying from to b and (b) cen- 
tered at {x, y) — (0, 6/2) with r varying from to 6/2. The ver- 
tical lines denote the threshold radius rt for tunneling within 
the first channel (see text). 



low energy: resonant transmission is larger from an outer 
lead to the inner lead which is closer to the opposite outer 
lead, where rotator modes can be accessed. 



C. Geometry-dependent mean transmission 

Having presented the zero-field spectral features and 
their origins, we proceed to investigate the overall multi- 
terminal transmittivity in dependence of the geometry of 
the setup by computing, for each set of parameters, the 
channel-integrated mean transmission coefficients 



in) 



(10) 



of the n-th transversal subband, which constitute a mea- 
sure of the overall response of the system upon an incom- 
ing wave in one of the leads. 

In order to confirm the role of rotator and librator 
modes in coupling symmetrically placed leads we first ex- 
plore the effect of a perturbing disk on the mean transmis- 



sion coefficients. In Fig. [3] the 



(n=l) 



Tij (the channel 



superscript is dropped for n — 1) are shown as a func- 
tion of varying radius of the disk, for three different ec- 
centricities. In Fig. |3] (a) the disk constitutes a circular 
recess of the elliptic boundary, and in Fig. [3] (b) it leaves 
the boundary of the billiard unperturbed, but partially 
blocks direct transport in its bulk. 



In case (a) the rotator modes are gradually destroyed 
with increasing r, because incoming waves from an outer 
lead are defiected on the concave part of the boundary 
into the billiard interior (an analogous situation is pre- 
sented in Ref. [2^ for a fixed rectangular cut). Similarly, 
the librator modes are destroyed since they rely on the 
focusing ability, predominantly around the y-axis, of the 
convex boundary. As a result, the direct pathways be- 
tween symmetric leads are depleted, leading to an abrupt 
decrease in T2i(r) and T^^ir) above a critical disk ra- 
dius, which is about Vc ~ w/3 for all e [note that the 
scaling r/b stretches the plots horizontally for larger e: 
rc(ei) « 0.10 6, rc(e2) « 0.08 6,rc(e3) « 0.06 6]. Be- 
low this critical disk radius, efficient guiding of rotator 
modes and focusing of librator modes can be considered 
robust to boundary pcrturbations4i For smaller e, 134 (r) 
remains substantial over larger r, since the larger bil- 
liard can support a larger number of direct scattering 
pathways between the inner leads. Deflection on the per- 
turbed boundary enhances scattering to asymmetrically 
placed leads, leading to increased T32 and Tia from zero 
for r ^ 0. 

In case (b) the librator modes are again rapidly de- 
stroyed by the disk, resulting in a decrease of T^^ir) sim- 
ilar to that in (a). Rotator modes sufficiently localized 
along the elliptic boundary survive up to some disk size 
and still connect the outer leads; especially in (b, £3), 
T21 (r) forms a characteristic high plateau until decreasing 
abruptly when the remaining free width between disk and 
boundary becomes smaller than the leadwidth w = a/5. 

The latter condition is met above a threshold radius 
rt given by (a) rt/6 = 1 — {5^/T~~e^)~^ and (b) 2rt/b = 
1- 



2(5-\/l — e^)"^, denoted by vertical lines in Fig. [3l at 
which the energy threshold for transport between oppo- 
site sides of the disk enters the first channel in the leads. 
When this threshold rises above k = 2 (for free width 
smaller than w/2), only tunneling through the constric- 
tions contributes to T21, T34 and 732, which then practi- 
cally vanish; in contrast, 713 is enhanced for large r, par- 
ticularly in case (b, ei). Furthermore, larger number of 
accessible eigenstates causes more resonant features in the 
transmission spectra, leading to increased fluctuations of 
the Tij in continuously varying geometry for larger bil- 
liards (row £3). 

Librator- and rotator-like eigenmodes of the billiard 
were shown to be necessary for high inner and outer lead 
transmission; nevertheless, their coupling to the leads fur- 
ther depends on the eccentricity of the unperturbed bil- 
liard for a given lead positioning. In Fig. 2] the Tij are 
plotted against the ratio b/a = y/1 — e'^ for different w. 
In order to access the limit of zero curvature (e = 1), the 
straight edge of the billiard is lowered by one leadwidth 
(sec dashed line in Fig.[T]). For the true semi-elliptic setup 
the features in Tij {b/a) are shifted to larger b (and thus 
smaller e) , so that the change in size is compensated and 
the corresponding eigenmodes remain approximately at 
the same energies. T21 and T34 overall increase with b, 
as rotator- and librator-like modes start to form which 
couple outer and inner leads; in contrast, 232 and T13 
overall decrease and possess a common broad minimum. 
For a certain eccentricity range a separation between in- 
ner and outer leads is thus possible, in the sense that 
cross-coupling between them (i.e. between an outer and 
an inner lead) is almost eliminated in zero magnetic field. 
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(a) Wo = a/3.5 



(b) W2 = a/7 




FIG. 4. (Color online) Tij for varying ratio b/a, with equidis- 
tantly attached leads of width (a) wo = a/3.5 and (b) W2 — 
a/7. The light (orange) dashed line shows the eccentricity e. 
The lower horizontal edge of the billiard here lies at y = —w 
(see Fig. [T]), and the upper edge is varied from straight {b — 0) 
to elliptic (0 < & < a) to circular {b = a), as illustrated by the 
inset pictures in (b). 



Further, depending on e, rotator-like modes leaking into 
the outer leads can interfere into a suppressed transmis- 
sion background, causing the characteristic wide dip in 
T21 for b « a/ 3. By decreasing w the transversal sub- 
bands are shifted up in energy and the wavelength of the 
incoming particle decreases relatively to the billiard size, 
so that a larger number of eigenstates is spanned within 
the first channel. This leads to increased fluctuations of 
the Tij in Fig. Ul^b), similar to those in Fig. [3] (row £3). 

Let us now investigate the overall transmission behav- 
ior in dependence of the placement of the inner leads, 
which determines their coupling to the different billiard 
eigenmodes. As the inner leads arc moved away from the 
center, there is a gradual crossover of the direct transport 
paths from librator- to rotator-type states. The ques- 
tion arises whether it is sufficient, or even necessary, to 
place the inner leads at the ellipse foci in order to achieve 
high transmission between them, as shown in Fig. [2j In 
the classical picture the separation of librators and ro- 
tators by the focal points is sharp, and in the limit of 
zero leadwidth, all trajectories coming in from one fo- 
cus are scattered directly (by only one reflection at the 
elliptic boundary) to the other, leading to unit transmis- 
sion. For a finite leadwidth, a portion of the incoming 
trajectories is scattered onto the straight segments be- 
tween the leads and eventually into an outer lead, so 
that the inner lead transmission is lowered from unity. 
In the quantum case, additionally, the spatial separation 
between librator and rotator modes is not sharp, espe- 
cially at the low energies considered; thus, however nar- 
row, the inner leads couple to both types of modes. Most 
importantly, though, the transmission coefficients highly 
depend on interference phenomena between the resonant 
states coupling to the leads. Even if the inner leads are 
placed close to the foci, where most eigenmodes possess a 
probability maximum,M multiple destructive interference 
between them may lead to low overall transmission. 

In Fig. [5] the variation of the mean transmission com- 
ponents with the inner lead displacement / is shown for 
different leadwidths and eccentricities. For large e, T34(Z) 
increases to maximum when the inner leads are next to 
the outer ones, though with transport dominated by libra- 
tor modes, since the foci lie within the outer lead openings 
{wi, ei). This trend is inverted for smaller e, where the 
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FIG. 5. (Color online) Tij in the unperturbed semi-elliptic 
billiard, for lead width (left to right column) wi = a/5, W2 = 
a/7, U13 = a/10.5 and eccentricities (top to bottom row) ei = 
0.78, 62 = 0.65, £3 = 0.35, as a function of the displacement I 
of the inner leads from the origin. The vertical lines show the 
position of the focal points //a = e for each setup. 



foci come closer to the origin and allow for the coupling 
to rotator modes. Then {wi, £3) T^i{l) is maximal for the 
inner leads close to the origin and decreases to minimum 
for large I, as a result of destructive interference. For nar- 
rower leads (columns W2 and ws) these features remain, 
with enhanced fluctuations like in Fig. 2] 

We indeed see that placing the inner leads close to the 
foci does not necessarily lead to high overall transmission 
between these leads, demonstrating the departure from 
the classically expected behavior of our billiard in the 
deep quantum regime. Some of the setups, in particular 
of Fig. [5] (u)3, £1), even exhibit a wide minimum in T34 
around I ^ f. 

^2i(0 overall increases with w and £, but remains 
largely unaffected by the variations in L It decreases 
slightly only when / is large enough for the inner leads 
to couple to the same modes as the outer leads, which 
causes a corresponding increase in T32 and Ti^. 

For the magnetic control of multitcrminal transmission, 
to be discussed in the following section, it is important 
to achieve high inner- and outer-lead zero-field transmis- 
sion, while cross-coupling is suppressed; we see that these 
conditions are met in our setup by combining small I and 
£ with relatively large w. 



D. Transmission in a magnetic field 

To understand how conductivity between terminals can 
be selectively manipulated with the magnetic field, we 
first analyze its impact on the transmission spectra and 
its interplay with the geometric properties. In the pres- 
ence of the field the resonant states accordingly shift in 
energy;^ while its influence on their phase modifies their 
coupling to the leads and the interference with other 
states. Therefore, the widths of sharp Fano resonances 
generally change, and the non-resonant pathways inter- 
fere into a different transmission background. If the 



field is very weak, the spatial distribution of the eigen- 
modes remains practically unaffected, as does their indi- 
vidual coupling to the lead openings. A drastic change 
in the overall transmission in a weak field can still take 
place, though, when a small number of leaking modes 
interfere4S For a stronger magnetic field the spatial dis- 
tribution of the states changes enough to generally yield a 
completely modified transmission spectrum. In the clas- 
sical picture the charged particle moving in the billiard is 
deflected into circular orbit of gyromagnetic (Larmor) ra- 
dius Tl = k/\B\, making the classification of trajectories 
into rotators and librators inapplicable. When the field 
strength is further increased, eventually becomes so 
small that the particle moves along skipping trajectories 
at the billiard edgesi^ The corresponding quantum scat- 
tering wave function is localized into edge states;^ which 
enable almost rcfiectionless transport along the boundary. 

In Fig. ini[a) the transmission spectra at (i) zero, (ii 
- iii) intermediate and (iv) high magnetic field strength 
are shown for the unperturbed billiard. As in Sec. IIII Bl 
the geometry is adjusted for high inner- and outer-lead 
zero-field transmission; in order to concentrate on mag- 
netically induced spectral changes, even wider leads are 
used, avoiding increased fluctuations from multiple inter- 
ference of billiard eigenmodes. The three first channels 
n = 1,2,3 are addressed, showing the typical stepwise 
increase of (maximal) transmission with n and the ef- 
fect of channel mode coupling. T21 and T34 are overall 
high in the first channel at S = [(a) (iii)], as previously 
discussed, but get modulated in the higher channels by 
increasingly wide dips and lowered background transmis- 
sion caused by multimode interference. The field strength 
for (a)(ii) is chosen to suppress inner- lead transmission, 
which enables the control of multiterminal transport, as 
will be shown in Sec. IIIIFl Waves coming in from lead 1 
are deflected onto the curved boundary, which decouples 
them further from the inner leads, reducing the width 
of the dips in T2i. is overall lowered, since waves 
coming in from lead 4 (containing the right ellipse focus) 
are no longer focused into lead 3: the librator modes are 
destroyed in the presence of the deflecting field. On the 
contrary, T32 = T14 is drastically increased and T12 (not 
shown) is accordingly reduced: the Larmor radius at this 
field is, approximately, the one needed to deflect classi- 
cal trajectories from lead 4 into lead 1 (or from lead 2 to 
lead 3) without reflection at the boundary, for energies in 
the flrst channel. This condition is not fulfilled anymore 
at the even higher field strength in (a) (iii), where edge 
states start to form. In the first channel, transmission is 
then favored to the next neighboring lead (T21, T34 and 
Tis = T42) and suppressed between other lead pairs. The 
complementarity between the multi-terminal coefficients 
is here clearly manifest in the coincidence of the dips in 
T34 and T13 = r42 with the peaks in r32 = T14. These 
resonances appear when the nodal pattern of interfering 
edge states (or a multiple of the diameter of the classi- 
cal skipping orbits) matches the distance of leads 2 and 
3 instead of adjacent leads (this behavior is described in 
detail in Ref. |331 for a 2-terminal billiard). In the second 
and third channel interference of a larger number of acces- 
sible modes enhances again fluctuations in the transmis- 
sion background. At very high held strength [(a)(iv)] the 
edge states lie so far apart in energy (and classically is 
so small) that plateaus of unit transmission to clockwise 
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FIG. 6. (Color online) Transmission coefficients T21 (red 
dasiied-dotted line), T34 (blue solid line), T32 (green dashed 
line), Tis (cyan dotted line), with offsets 3, 2, 1, respectively, 
as a function of k within the first three transversal subbands 
n — 1,2,3 of leads of width w = wo — a/3.5, attached 
equidistantly with I = {a — w/2)/3 = 0.84/, (a) for the 
unperturbed semi-ellipse billiard with e = 0.35 and (b) for the 
same geometry with a disk of radius r = b/2 — 2w/3 centered 
at (0,6/2), at magnetic field (i) B = 0, (ii) B = 0.002, (in) 
B = 0.005, (iv) B = 0.010, with B = Bz pointing outwards 
from the billiard plane. 



subsequent leads appear for energies in the lowest magne- 
toelectric (Landau) subband. When more edge states are 
energetically accessible we observe a difference between 
scattering at smooth and sharp lead openings. Diffrac- 
tion at the sharp edges causes mixing and interference of 
the different edge states,— leading to oscillations in T34, 
T32 and T13. Only the T21 coefficient exhibits perfect 
transmission even at higher energy, since the edge states 
adiabatically follow the smooth elliptic boundary from 
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FIG. 7. (Color online) 
Mean transmission co- 
efficients Tij as a func- 
tion of the applied mag- 
netic field B for the bil- 
liard of Fig. [S] with in- 
creasing disk radius r 
in the n — 1 subband 
(top to bottom panels 
in left column) ; for zero 
and maximal r (top 
and bottom row, re- 
spectively) also the n = 
2 and n = 3 subbands 
are shown. Each plot is 
of height max[fjf] = 
n. The light (orange) 
dashed line in the bot- 
tom row panels shows 
the scaled Larmor ra- 
dius tl/ci at the n-th 
channel's center k = 
n + 1/2. 



lead 1 to 2. The stepwise increase of transmission with 
n, most pronounced in T21, is shifted to higher k with 
increasing field strength, following the threshold energies 
of the magnetic subbands. 

As in Sec. IIII C[ a disk inside the semi-elliptic billiard 
drastically changes its transmission properties by block- 
ing direct transport paths between the leads. In Fig. [Bib) 
the field dependent spectra are shown for this setup, with 
the disk leaving constrictions of minimal width 2w/3 with 
the boundary. The transmission threshold is then essen- 
tially shifted from n to 3n/2, as can be seen by replacing 
w with 2w/3 in Eq. (0). Thus, T21, T34, or T32 > n - 1 
below K = 3n/2 results from tunneling of the wave func- 
tion through the constrictions. Distinct resonant tunnel- 
ing peaks are seen below this threshold for n = 2 (that 
is, 1.5 < K < 3) in and T34 at high field strength 
[Fig. IHb)(iii) and (iv)], mediated by edge states that are 
localized on the disk edges (similar to the states leading to 
sharp reflection resonances, due to different geometrical 
setup, in Ref. l49l ). The effective resonator length is, ap- 
proximately, the mean periphery 2n{r + w/3), leading to 
the observed peak spacing Ak w 0.247 in (b)(iv). On the 
other hand, scattering upon the disk favors transmission 
between leads on the same side of it, so that T13 overall 
increases and qualitatively approaches the unperturbed 
case for strong fields. 

In Fig. [7] the four independent channel mean compo- 
nents Tij are plotted as a function of B for different radii r 

— (2) 

of the inserted disk [rows (i) to (v)]. The T>j and T^j arc 
shown for r = (unperturbed billiard) and r = 6/2 — w/3 
(almost divided billiard), where also the Larmor radius 
rL{B) = -k/w ■ {n + ^)/\B\ at each channel center is plot- 
ted to show the field impact on the classical trajectories. 
In absence of the disk [row (i)], T21 is close to unity for 
i? ^ 0, falling off slowly at large B as the magnetic thresh- 



old enters the channel [see Fig. [nija)(iv)]. For i? < it 
decreases abruptly with field strength when < a, since 
the incoming waves in lead 1 are deflected away from the 
elliptic boundary. 734, which is also large at B = 0, de- 
creases to a prominent local minimum at B « +0.002, 
corresponding to the spectrum in Fig. [ni[a)(ii). At the 
minimum, a large portion of the wave coming in from 
lead 4 is deflected into lead 1, leading to a corresponding 
maximum in T'14 [with Ti4^{B) = T'4i(— B) = T^2{B), see 
Table H], which is otherwise close to zero. It is this dras- 
tic change at the intermediate field strength B w +0.002, 
different for each mean transmission component, that will 
serve as a key property to enable multiterminal transport 
control for the geometrical parameters used. At even 
stronger field the waves from lead 4 follow edge states 
directly into lead 3, so that T'34 increases again, while 
Ti4 = T^2 decreases. For i? < 0, T34 remains high over 
the local minimum of T21 around B = —0.002. before it 
too falls off for stronger fields. Tis, on the other hand, 
remains close to zero for all i? < 0, increases with the 
field strength for i? > 0, and finally follows T'34 in the 
edge state regime, where the pathways 4 — >■ 3 and 3 — > 1 
are almost equivalent along the boundary: indeed, the 
r34 and ri3 spectra practically coincide in Fig. [BJa)(iv). 

On all Tij(B) curves, though more visible in T21 and 
T'34, relatively small fluctuations in B are superimposed, 
which can be regarded as generalized collective AB oscil- 
lations from interference between spatially extended leak- 
ing states: the oscillations in Tij{B) at each k, add up to a 
large-scale oscillation of the channel average. The charac- 
teristics of Tij remain qualitatively the same in the higher 
channels (with maximum = n), mapped onto a larger B- 
scale: at higher energy larger fleld strength yields the 
same Larmor radius and similar variations as in Tij(B). 
From the above we see that, depending on B, overall 
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transmission is favored from each lead to certain other 
leads and suppressed to the rest. We will address this 
possibility for directed multiterminal transport in detail 
in Sec. HITFI 

The modification of the Tij{B) profiles by the perturb- 
ing disk is shown in Fig. ^\) to (v), where its radius r 
is increased so that transmission between leads on oppo- 
site sides is suppressed, as previously described [see Figs. 

andUb)]. Thus T2i(B), Tm{B) and ^32(5), al- 
though retaining their trends, gradually decrease to zero 
for any B when the constrictions become narrower than 
the leads. In contrast, T13 increases with r at i? = 0, 
as seen also in Fig. Eljb), and remains large at strong 
B > Q. Also at strong B < 0, though, Tia increases 
with r, because the edge states (now clockwise deflected 
classical orbits) can guide the particle from lead 3 onto 
the disk edge and then onto the elliptic boundary to the 
left of the disk, which it follows into lead 1. Interest- 
ingly, for large enough disk [as in Fig. where the 
constriction width is if/ 3] the 4-terminal billiard is ef- 
fectively divided into two 2-terminal billiards for n = 1, 
so that transmission between leads on the same side of 
the disk becomes symmetric in B: Ti^{B) = Ti3(— _B) 
and T2i{B) = T2i{—B), as a consequence of the sum rule 
Eq. © and the symmetry relation Eq. ([5]) for the spe- 
cial case of a system with two leads This 2-terminal 
symmetry is not present in the higher channels [n = 2, 3 
in row (v)], where the smaller transversal wavelength en- 
ables transport through the constrictions. 

Conclusively, the disk reduces the difi^erence between 
the (independent) Tij{B), thereby weakening the control- 
lability of mutliterminal transmission. Nevertheless, an 
appropriate disk-like blocking potential switches output 
from lead 2 to lead 3 (from lead 3 to lead 2) with input in 
lead 1 (in lead 4) at strong i? > 0, which, in the context 
of directed transport, constitutes an additional (electric) 
switching mechanism based on the geometry-independent 
behavior of edge states. 



E. Bent coupled wires 

Magnetically induced directed transport in our setup 
was shown to rely on the convexity of its boundary, which 
enables highly transmittivc pathways at low field strength 
by the formation of librator- and rotator-like modes. 
These coexist in the semi-ellipse, so in order to separately 
investigate the role of boundary-localized modes in the 
overall transmission, and the impact of their interference 
with gradually upcoming bulk modes, we consider the fol- 
lowing setup: two circularly bent parallel quantum wires 
coupled through a smooth opening of adjustable width 
d at the bend [sec sketch in Fig. [8l[a)]. The Tij{B) arc 
plotted in Fig. |Sl^a) without the opening (top) and for in- 
creasing opening width (second from top to bottom), and 
are compared to the case of straight coupled wires^^iS in 
Fig. IHl^b), where curvature is absent and additional y- 
symmetry is present. With no opening the transmission 
in the two disconnected bent wires at i? = is almost per- 
fect: T21 and T34 depart from unity only due to narrow 
resonances caused by the curvature of the wires, which ef- 
fectively induces an attractive potentiali^ is slightly 
smaller than 734, because the resonances for the longer 
bent part lie closer in k. Both slowly decrease at stronger 
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FIG. 8. (Color online) Tij(B) for two parallel wires at dis- 
tance equal to their common width w — a/3.5, coupled by a 
smooth opening of width d and circular edges of radius w/2, 

(a) with the wires bent by an angle tt across the coupling and 

(b) in straight configuration (parallel and equidistant to the 
a:-axis). In the top plots there is no opening, and then the 
opening width is increased in steps of ui/4 from d = to/4 to 
d = 3w (second from top to bottom, with transmission offsets 
decreasing by one), the latter yielding in (a) the semi-circular 
billiard with smooth lead openings. 



fields, as the magnetoelectric subband threshold rises [like 
in Fig.[nra)(iv)]. 

When the wires couple, the 2-terminal S-symmetry of 
the isolated wires is broken. As the opening is widened, 
the rotating modes connecting the outer and inner leads 
increasingly interfere with states that extend into the 
opening, leading to enhanced back-scattering dips in the 
corresponding transmission coefficients and a following 
decrease in T21 and T^i. The role of rotator- and librator- 
like modes is here manifest in the inner-lead transmis- 
sion as d is varied at _B « 0. T34 is high for d « 0, 
where transport is dominated by rotator-like modes, de- 
creases at intermediate d, where the rotators are de- 
stroyed by the opening in the inner convex boundary, and 
increases again for large d (> 6w/A), where the outer con- 
vex boundary focuses the wave function into librator-like 
modes. At appropriate strength, the magnetic field favors 
transport between terminals on either side of opening by 
deflection of the particle orbits. T32 increases to a local 
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FIG. 9. Multiterminal conductance coefficients gij at scaled Fermi energy kf ~ w/n ■ \/2Ef — 1.425 and temperature = 
100 mK, for B — 0, B = ±0.002 and B — ±0.010. In each 4x4 block, row i labels the input and column j the output current 
terminal. 



maximum at intermediate B > 0, which rises for larger 
d. Ti3, on the other hand, rises significantly only at high 
B > where the particle is guided by edge states. 

For maximal opening d = 3w the setup becomes the 
4-terminal semi-circular billiard, now with smooth lead 
openings. Controllability of output terminal, as described 
in Section Fill Dl then becomes optimal, and the Tij{B) 
profiles are very similar to the ones in Fig. ^\). This 
shows that the smoothness of the lead openings, as well 
as the small change in eccentricity and size, although 
clearly affecting the dynamics in the scattering system 
and thereby the detailed spectral features, leave the over- 
all field dependence qualitatively unchanged. 

In Fig. |Hl[b) the strong field asymptotics of the Tij{B) 
for the straight coupled wires coincide with those of the 
bent wires, since transport through edge states is rather 
affected by the topology, and not by the geometry, of the 
scatterer. The first obvious difference here is that the 
additional reflection symmetry about the x-axis renders 
the components T21 and 734, whose difference was central 
in the discussion so far, identical. That is, regardless of 
the field strength, high transmission between the outer 
leads can never be combined with low transmission be- 
tween the inner leads, as it can for the bent wires. Also 
due to spatial symmetry, T32 must now be symmetric 
in B, and its broad peak around B = ±0.002 for the 
bent wires is shifted to B = Q. As a result, crossed- lead 
transmission (between leads 2,3 or 1,4) can no longer 
be switched from high to low by inverting the field. At 
weak fields, T^4^{B) varies similarly for bent and straight 
wires for small d < 6^/4; for large d though, the straight 
wires, unlike the bent wires, yield low 234 = T21 around 
i? = 0, because of the absence of modes that are focused 
at the lead openings or guide the incoming wave along the 
boundary. We conclude that the increased symmetry of 
the straight leads and the absence of the convex boundary 
reduce the possible combinations of magnetically induced 
transport directions between the terminals. 



F. Directed conductance 

In the discussion so far we have utilized the channel- 
integrated mean transmission [Eq. ([TU|) ] as a tool to 
compactly describe the average response of the trans- 
mission to parameter changes. The actual measurable 
conductance coefficients gij connecting the current flow- 
ing inwards at terminal i with the voltage differences 
to all other terminals j are, in the linear response 
regime at temperature 0, given by the Landauer-Biittiker 



formula;^ in units of the (spin-degenerate) conductance 
quantum 2e^/ft., 

/+00 
T,,{E)F{Q,EF\E)dE (11) 
-00 



where 



1 ^(E-Ef, 

sccli [ _ „ 1 (12) 



with f{EF\E) being the Fermi distribution for given 
Fermi energy Ef- Thus, Qij essentially equals the ther- 
mally averaged transmission around the electron Fermi 
energy, with a width proportional to 0, and coincides 
with Tij at e = 0. 

The output controllability described in Sec. IIII Dl can 
be enhanced in terms of the conductance at low temper- 
ature and low Fermi energy, where transmission features 
are resolved in a smaller K-range than the whole chan- 
nel. While the local maxima (minima) of 232 and 234 in 
Fig.[71[i) suggest a maximal efficiency for magnetically di- 
rected transport of about 50% in the unperturbed elliptic 
setup, adjusting temperature and Fermi energy appropri- 
ately can yield corresponding conductance maxima (min- 
ima) close to unity (zero) at these field strengths. To pro- 
vide an example, we choose the Fermi energy correspond- 
ing to K ^ kf ^ \fTEF • w/tt = 1.425 [cf. transmission 
spectrum in Fig. IHKa)] and a temperature = 100 mK. 
The resulting conductance coefficients for the above de- 
termined optimal field strengths are shown in Fig. 15] as 
grayscale cells for the individual gij ordered like in Ta- 
ble m The diagonal elements gu following from Eq. ([TT|) 
do not contribute to calculated currents, but indicate the 
degree of reflectance for ballistic transport and show the 
depart from unity of the sum of conductances from or 
to the other (^ i) terminals. For the chosen parame- 
ters the conductance coefficients practically reach unity 
(black cells) for specific terminal combinations i -> j and 
practically vanish for the rest (white cells), depending on 
the direction and strength of the field. In particular, for 
each input lead i, the output can be switched selectively 
to any lead j ^ i hy appropriately tuning B. This relies 
on the above discussed interplay of geometry and mag- 
netic field effects. For large \B\ (see B ~ ±0.010 blocks in 
Fig. [9]) edge states form and conductance is determined 
by the topology of the boundary (in general directed edge 
state transport can also be implemented, but with finite 
potential barrier a^"^i^^'^^ ). At zero and intermediate \B\ 
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the output is governed by interference of spatially ex- 
tended scattering states leading, for this billiard, to over- 
all high transmittivity between certain non-neighboring 
terminals, as seen in the B = and B = ±0.002 blocks 
in Fig. M 

This shows that the current can be efSciently directed 
from a given input to a selected output terminal. Espe- 
cially for zero and intermediate magnetic field strength it 
is even possible to construct a controllable cross-junction 
for the currents: choosing, e.g., leads 1 and 3 as input 
terminals and leads 2 and 4 as output, current is flowing 
from lead 1 to 2 and from lead 3 to 4 for i? = 0. For 
B = -1-0.002 current is flowing mostly from lead 1 to lead 
4 and from lead 3 to lead 2, thus exchanging the directed 
connections between the terminals with applied magnetic 
field. 

For higher Q the gij{Q;Ep) in Fig. IH] arc generally 
shifted to more intermediate values, because F{Q, Ep; E) 
is broadened and additional spectral features arc included 
in the integration around Ep in Eq. ([TT]) . On the other 
hand, for very low the transmission spectrum is highly 
resolved by F(9, Ep; E) and the gij become sensitive to 
small changes in Ep. These dependencies are not present 

~(n) 

in the channel means , which therefore serve to esti- 
mate the field strength values suitable for directed trans- 
port in a given billiard; with adjusted temperature and 
Fermi energy, the effect can then be optimized by slightly 
modifying the transmission features through fine-tuning 
of B, in order to controUably obtain maximal and mini- 
mal gij. 

IV. SUMMARY 

We have investigated the ballistic transport proper- 
ties and low-temperature magnetoconductance of a 4- 
terminal semi-elliptic quantum billiard. Analyzing the 
strong dependency of the transmission on the magnetic 
field and geometrical parameters like eccentricity of the 
semi-ellipse and placement and width of the leads, we 
have shown how electrons can effectively and controllably 
be guided from one input lead to any other output lead, 
including the cross-switching of output from combined 
input. 

At zero field strength electrons are guided between the 



two outer or the two inner leads by modes corresponding 
to classical rotator or librator trajectories, respectively. 
The role of rotator and librator modes has been clari- 
fied by introducing a perturbing disk, thereby destroying 
modes selectively. In this context we have shown that 
the transmission is robust under small perturbations and 
that the disk offers additional possibilities of transmis- 
sion control. The efficiency of the selective transmission 
has been optimized with respect to the geometric param- 
eters. Here it turns out that the optimal position of the 
inner leads deviates from the classically expected one at 
the focal points of the ellipse and that an eccentricity of 
about 0.35 for equidistant terminals is optimal. 

The application of the magnetic field allows to control- 
lably change the conductance coefficients: at appropriate 
field strengths electrons coming from one input lead can 
be guided to any other output lead. This results from 
the deflecting effect of the field at intermediate strength. 
Further, by increasing the field strength to the edge state 
regime, conductance is mostly determined by the topol- 
ogy of the billiard. In this regime sharp edges cause mix- 
ing and interference of multiple edge states, leading to 
oscillations in the transmission spectra, while a smooth 
boundary guides the individual edge states without mix- 
ing. 

To further examine the role of the elliptic boundary, i.e. 
of the existence of rotator and librator modes, we have 
investigated the transmission through a pair of bent cou- 
pled quantum wires and compared it to the topological 
equivalent setup of two straight coupled wires. The bent 
wires show a degree of control superior to the straight 
setup, but the semi-ellipse allows for the highest degree of 
conductance control. While the channel-averaged trans- 
mission clearly shows the possibility of directed transport, 
very high switching efficiency is achieved at low but re- 
alistic temperatures and appropriate values of the Fermi 
energy. The semi-ellipse qualifies as a magnetically con- 
trollable cross-junction for ballistic quantum transport. 
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